Tip induced doping effects in local tunnel spectra of graphene 
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We report on tip induced doping in local tunnel spectra of single layer graphene (SLG) with tun- 
able back-gate using room temperature scanning tunneling microscopy and spectroscopy (STM/S). 
The SLG samples, prepared on silicon dioxide surface by exfoliation method and verified by Raman 
spectra, show atomically resolved honeycomb lattice. Local tunnel spectra show two minima with a 
clear evolution in the position of both with doping by the back gate. A similar variation in spectra is 
also observed spatially due to charge inhomogeneity. With doping the two minima move in opposite 
directions with one showing nearly a square root dependence and the other a linear dependence 
on gate voltage. We explain these features as arising from the STM tip induced and bias voltage 
dependent doping in SLG. 



The electrical transport in graphene and its depen- 
dence on carrier concentration has revealed some spec- 
tacular physics in this Dirac-like two dimensional elec- 
tronic system The carrier concentration in single 
layer graphene (SLG) has mainly been controlled by a 
silicon back gate separated from graphene by a 300nm 
thick Si02 layer. Two gates, one on top in addition 
to the bottom gate, have been used to create a tunable 
gap in the bilayer graphene [2] . Recently, effects of local 
doping on bulk transport have also been reported us- 
ing scanning gate microscopy [3] and scanning tunneling 
microscopy (STM) [4] which use the microscope tip to 
induce doping. Scanning tunneling microscopy and spec- 
troscopy (STM/S) is an ideal tool for locally probing the 
tip induced doping. Here we report on tip induced dop- 
ing effects in local tunnel spectra on SLG prepared by 
mechanical exfoliation on Si02/Si surface. Local tunnel 
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FIG. 1: a. Optical image of few layers graphene (FLG) with a 
portion (marked as 'Graphene') of SLG on Si02 b. The FLG 
flake after deposition of Cr/Au contact leads, c. Raman 
Spectra of SLG region in ambient conditions using laser light 
of wavelength 514 nm, spot size 1/im and power 0.64 mW. The 
inset shows the details of the 2D peak centered at 2681 cm -1 , 
d. STM measurement schematic of this graphene device. 



spectra exhibit two minima with the two moving in op- 
posite directions along the bias axis with increasing gate 
voltage. We explain the two minima and their evolution 
with gate voltage in terms of the bias voltage dependent 
doping by the STM tip. 
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FIG. 2: a. Topographic image of SLG (560x560 nm 2 at 
0.4V and O.lnA). The inset shows topographic profile along 
the white line shown in the image b. Atomically resolved 
constant-current STM image (2.4x2.1 nm 2 at 0.4 V) of SLG. 
c. Topographic image of SLG (161x161 nm 2 at 0.4V and 
O.lnA). d. Tunneling conductance at various positions shown 
in c. 

The SLG samples were prepared on 300 nm thick insu- 
lating Si02 film on doped silicon substrate by mechanical 
exfoliation of graphite flakes [5]. SLG as a part of few lay- 
ers graphene (FLG), size(^ 50/im x 18/im) was optically 
identified [6] as shown in FigQJi. For connecting leads we 
deposited ~1 mm wide and few mm long strips of lOnm 
Cr topped by 100 nm Au by masking larger portion of 
graphene using a tungsten wire of 25 fim diameter (see 
FigHb)- Such mechanical masking [?J minimizes the con- 
tamination from the liquid chemicals used in lithographic 
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techniques. The /i- Raman spectrum in FigJT] shows the 
G peak at 1580 cm -1 , 2D peak at 2681 cm -1 and the 
minor D peak at 1339 cm -1 . This confirms this region 
to be SLG 0. 

The STM data reported here were taken at room tem- 
perature using a homemade STM in a cryo-pumped vac- 
uum chamber with 10 -4 mbar pressure. The gate voltage 
(V g ) was applied to the silicon substrate with 270 kQ se- 
ries resistance. For STM measurements, as depicted in 
the schematic in Fig. [T]i, we apply the bias voltage (Vb) 
to the sample while the tip stays at (virtual) ground po- 
tential. A homemade 2D nano positioner [9[ was used 
for aligning SLG with the STM tip under an optical mi- 
croscope. Electrochemically etched and HF treated [loj 
tungsten wire of 0.25 mm diameter was used as the STM 
tip. The apex radius of the tip was found to be in 20- 
50nm range from electron microscopy. The tunnel spec- 
tra were acquired using an ac modulation of amplitude 
20 mV with the dc bias voltage. The junction resistance 
was kept as 4GS1 (400mV bias and lOOpA tunnel current) 
for all the local spectra. We have checked reproducibility 
of the spectra with several tungsten tips on different SLG 
samples but the data presented here are from the sample 
shown in FigJTJ 

FigE^i shows the topographic image of a large area 
(560x560 nm 2 ) of the SLG region of the sample. The yel- 
low and dark regions are bumps and troughs either due to 
charge inhomogeneity or ripples as observed in previous 
STM studies [TT], [l2| . Inset of this image also shows the 
topographic profile, with a roughness of ~0.2nm, along 
the white line marked in the image of Figj2^i. Figj2]3 
is the atomic resolution image showing the honeycomb 
lattice structure of SLG. Fig|2t shows a topographic im- 
age of area 161x161 nm 2 with the tunnel spectra at the 
marked locations shown in the Figj2]i. These spectra are 
taken at zero gate voltage (Vg). All of the spectra show 
two minima contrary to the expected 'V shape DOS for 
SLG. The first minima (closer to zero bias) has less curva- 
ture than the second minima. The position of the second 
minima on the bias axis varies between 0.4 and 0.5 V 
at different locations. The location of the two minima 
varies as one moves to different locations of the sample 
due to local charge density [12| as we discuss later. 

We can understand the two minima if we incorporate 
the doping arising from the tip bias voltage as follows. 
Using the linear dispersion for graphene near Dirac point 
[l|, i.e. ek± = =L/2/up|k| with vp = 10 6 m/s as the Fermi 
velocity, and the 2-fold spin and 2-fold valley degeneracy, 
the DOS of SLG is given by, 
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Here, Ep is the Fermi energy measured from the Dirac 
point, and determined by the carrier density, and E is 
the energy measured with respect to Ep. For a given 
number densi ty of e lectrons, n, and using Eq. [I] we find 
Ep = Tivpn^/it /\n\. Assuming a parallel plate capacitor 
configuration (see Fig. [T]i) of SLG with both the bottom 
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FIG. 3: a. Theoretical tunneling conductance at T — OK for 
different gate voltages. The inset shows a comparison between 
T = OK and 300K for V g = 20V. b. Variation of primary and 
secondary minima with V g . The bottom inset shows the low 
bias region of the spectra at small V g values showing the fine 
spectral features. The top inset shows the minima locations 
for small V g and Vb values. All the plots are for /3 — 75 and 
7 = 0.032V 1/2 . 



gate electrode and the STM tip [l3[, we can write, 
Ke 0i 



ed 



<v g - m). 



(2) 



Here, /3 = with d as the oxide layer thickness, z is the 
tip-sample separation and K (« 4) is the dielectric con- 
stant of Si02- e is the magnitude of electronic charge, eo 
is permeability of free space, V g is the gate voltage and Vb 
is the bias voltage applied to the sample keeping the STM 
tip at zero potential. From here we see that although the 
magnitude of Vb is smaller, by about two orders of mag- 
nitude, than V g , the factor j3 (= 75 for z = lnm) makes 
the two terms comparable to each other. Thus the tip 
induced doping effects should be easily observable. 

Using Eq. [2] for n to find Ep, which in turn 
can be used in Eq. [T] to get DOS as, N(E) = 
-pv b 



E + 7e- 



Here, 7 
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0.032 V 1 / 2 for v F = 10 6 m/s, d = 300nra and K = 4. 
Thus the tunnel conductance at zero temperature is given 

by 



G(V b ,T = 0) = ea 



v\v 9 -m\ 



(3) 



Here is a is a constant. Since the tunnel conduc- 
tance at finite temperatures is proportional to the con- 
volution of DOS and negative of the derivative of the 
Fermi function, f(E,T), [13, we get, G(V h ,T) = 



aj_ 



E 



1 e(V g -pV h ) 



y/\Vg-PVb\ 



df(E-eV b ,T) 

dE 



) dE. Fig. 



shows the calculated spectra for /5 = 75 and 7 = 0.032 
V 1 / 2 with a minima other than the Dirac point and with 
a clear evolution with V g . At room temperature the Dirac 
point also appears as a minima due to thermal smearing 
and local disorder. These two minima are consistent with 
the experimental spectra. 
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FIG. 4: Experimental spectra at different V g values taken 
at a particular location on SLG. The spectra have been off- 
set vertically uniformly for clarity, b. Behavior of the two 
minima with V g . Scattered data show experimental observa- 
tion. Solid and dashed lines show theoretical expectation with 
P — 75 and 7 = 0.032V 1 / 2 . The inset shows three spectra (at 
V g =+30, and -30 V) without offset. 

Let us analyze the tunnel spectra features at T = 0, 
given by Eq. [3l Finite temperature only smears the spec- 
tra retaining the essential features. Due to square-root 
dependence of Ep on n, Ep has a sharp jump with n 
near n = 0, i.e. as Ep passes through the Dirac point. 
This happens at V g = fiVb and gives rise to a minima in 
G(T4,0) at \Vb\ close to but slightly larger than \V g /P\. 
The exact location of minimum at T = is found to be 

yrnin = Vg_ + l_f^_ for > £jl = L27 y for q = 75 and 

7 = 0.032 V 1 / 2 . Other than these minima the Dirac point 
exists at V b = ~^\(h 2 + 7\/ 4 l^l + ^V)- There is 
more structure in the conductance spectra for \V g \ < ^-j- 

and at low bias (|T4| < = 38 mV); however, it is not 
possible to observe these low bias features at room tem- 
perature due to limited energy resolution. Moreover, the 
low bias region may be significantly affected by disorder 
[Hj and electron-phonon interactions [16]. The spectra 
on our sample do not show a sharp cusp corresponding 
to Dirac point due to significant disorder as seen from 
the STS studies and also due to finite temperature. It is 
well known that SLG on Si02 has significant charge in- 
homogeneity due to local dopants that shifts the Ep [121 
away from the Dirac point. In addition local dopants can 
also give rise to changes in DOS [17j. In general, some 
feature in the spectra should be observable whenever Ep 
passes through a sharp feature in DOS. 

We also note that the second minima, when far from 
the Dirac minima, does not get thermally smeared by 
temperature as seen in the inset of Figj3^i. This happens 
as there is no actual minima in the DOS at this energy. 
This minimum comes from the linear dispersion at Dirac 
point and will occur only for SLG and thus the second 
minima will not be present in spectra taken on more than 
one layer of graphene though the location of the first 
minima may be affected by the tip induced doping. 

FigH^i shows the dependence of the spectra on gate 
controlled doping. These spectra were acquired at the 
same point of an image (not shown). We see two min- 



ima in these spectra with a clear evolution with the gate 
voltage. With increase of carrier concentration of ei- 
ther type the two minima move in opposite directions 
as expected. The exact locations of the minima were 
obtained by fitting a quadratic function in a small bias 
range (^100mV) around the minima and are indicated 
by the arrows. Fig|4)3 shows the comparison between the 
observed location of the two minima and the calculation 
with z = lnm. We see a good agreement between exper- 
iment and theory. Slight deviation of the minima from 
the expected behavior may arise from changes in spectra 
due to disorder [l5| and impurities [l7| . 
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FIG. 5: a. a. Spectra at different V g values (10V interval) 
taken at another location of SLG. The spectra have been off- 
set vertically uniformly for clarity, b. Behavior of the two 
minima with V g with the scattered data showing experimen- 
tal observations and solid line showing the calculated minia 
with z = 0.55nm and vf — 10 6 m/s. The dashed line shows 
the calculated minima z — 0.55nm and vf = 5 x 10 5 m/s. The 
inset shows two spectra from a at gate voltages 10 and -60 V. 

As discussed earlier, the two minima have significant 
spatial inhomogeneity arising from variations in local 
doping. Fig. [5] shows the V g dependence of the two 
minima at a different location showing this region to 
be electron doped with a carrier density equivalent to 
V g ~ —25V. In this case we obtain a good agreement 
with theory for the second minima using a tip-sample 
separation of 0.55 nm and incorporate the effect of elec- 
tron doping due to charge impurity. The observed shift 
in primary minima is smaller as compared to the theory. 
We can fit both the minima, as shown in Fig. [5j by using 
a 50% reduction in vp which may arise from the disor- 
der [l8[ as this SLG portion is close to dopants causing 
disorder. 

In conclusion, we have done a STM study of gated 
single layer graphene films prepared by mechanical exfo- 
liation. The single layer graphene shows a topography of 
0.3 nm over the 500 nm length scale. STS measurement 
at room temperature shows two minima whose shape and 
position changes with doping. We explain these features 
in detail on the basis of tip induced doping effects. 
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